Non-universal behavior of helicity modulus in a dense defect system 
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Extensive Monte Carlo simulation has been performed on a 2D modified XY model which behaves 
like a dense defect system. Topological defects are shown to introduce disorders in the system 
which makes the helicity modulus jump non-universal. The results corroborate the experimental 
observation of non-universal jump of the superconducting density in high-T c superconducting films. 
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The two-dimensional (2D) XY model of planar spins 
has been the subject of intense interest since the last four 
decades. The 2D XY model was believed to be with- 
out a phase transition for a long time until Kosterlitz 
and Thouless (KT) proved that a phase transition in- 
deed occurs and clarified its topological nature [l|, 
KT transition is a continuous phase transition of infi- 
nite order from a high-temperature isotropic phase with 
exponential decay of spin-spin correlation functions to a 
low-temperature pseudo-long-range or quasi-long- range 
order (QLRO) phase with power law decay of spin-spin 
correlations. KT pictured this transition in terms of vor- 
tex unbinding. In the low-temperature phase the charges 
(vortices) are bound together into dipole-pairs while in 
the high-temperature phase some dipole-pairs are bro- 
ken. KT theory Q leads to the famous universal jump 
prediction 0] of the helicity modulus [H or equivalently 
of the superfluid density @, Q . Kosterlitz RG equations 
for the 2D Coulomb gas (CG) are constructed in the low- 
temperature phase and are valid in the limit of small par- 
ticle densities [6| . Later Minnhagen has suggested on the 
basis of a new set of RG equations for the 2D CG that 
the conclusions based on Kosterlitz RG equations may 
break down for larger dipole-pair fugacities @, As 
a result charge unbinding transition with non-universal 
jumps may, in principle, be possible and it was shown 
that for higher particle densities, the charge-unbinding 
transition is first order Later Zhang et.al on the 
basis of sine-Gordon (SG) field theory, showed that the 
nature of the transition in the dense 2D classical CG is of 
discontinuous first order type. The evidence of a first or- 
der phase transition for higher particle densities was sup- 
ported by Monte Carlo (MC) simulations [llrEl- Such 
first order phase transitions are related to high-T c super- 
conductivity [l4j |. Leemann et.al, in their experimental 
measurements of the inverse magnetic penetration depth 
A -1 in thin films of YBa 2 Cu 3 07 [lfij], observed that the 
jump of the superconducting density (or equivalently the 
helicity modulus) did not obey the universal prediction of 
KT theory. These systems thus cannot be described by 
conventional XY model. Mila attempted to under- 



stand the the non-universal jump of the superconducting 
density in thin films of high-T c superconductors in terms 
of XY model with a modified form of interaction poten- 
tial. Such kind of model was first introduced by Domany 
and co-workers [ijj ■ They introduced an extension of the 
2D XY model where the classical spins (of unit length), 
located at the sites of a square lattice and free to rotate 
in a plane, say the XY plane (having no Z-component) 
interact with nearest-neighbors through a modified po- 
tential 
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where is the angle between the nearest neighbor spins, 
J is the coupling constant and p 2 is a parameter used 
to alter the shape of the potential, or in other words, p 2 
controls the nonlinearity of the potential well, although 
variation in p 2 does not disturb the essential symmetry 
of the Hamiltonian. For p 2 = 1, the potential reproduces 
the conventional XY model while for large values of p 2 
(say p 2 = 50), the model behaves like a dense defect sys- 
tem [l8| and gives rise to a first order phase transition as 
all the finite size scaling (FSS) rules for a first order phase 
transition were seen to be nicely obeyed 19]. The first 
order phase transition is associated with a sharp jump in 
the average defect pair density pH, Ho| . The change in 
the nature of the phase transition with the additional pa- 
rameter p 2 is in contradiction with the prediction of RG 
theory according to which systems in the same univer- 
sal class (having same symmetry of the order parameter 
and same lattice dimensionality) should exhibit the same 
type of phase transition with identical values of critical 
exponents. In this context, I refer to the work of Curty 
and Beck [2l[ who showed that in three dimension (3D), 
continuous phase transition can be preempted by a first 
order one. 

Enter and Shlosman finally provided a rigorous proof 
[13, Hl| of a first order phase transition in various SO(n)- 
invariant n- vector models which have a deep and narrow 
potential well. The model defined by Eqn. (JlJ is a mem- 
ber of this general class of systems. Moreover, Enter and 
Shlosman argued that in spite of the order parameter in 
the 2D systems with continuous energy spectrum being 
predicted to vanish by Mermin- Wagner theorem [24[ , the 
first order transition is manifested by the long range or- 
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der in hi gher -order correlation functions. Recently Sinha 
and Roy [19j verified this argument by numerical simula- 
tions and showed that while the lowest-order correlation 
function decays to zero, the next higher-order correlation 
function has a finite plateau. Later they investigated the 
role of topological defects on the phase transition exhib- 
ited by the model described by Eqn. ([1]) by means of 
extensive MC simulations and observed that the system 
appears to remain ordered at all temperatures when con- 
figurations containing topological defects are not allowed 
to occur [l8j . 

However, the connection of spin-stiffness (helicity mod- 
ulus) with topological defects has not yet been studied 
for systems exhibiting first order transition. The present 
Communication aims at studying this aspect of phase 
transition which allows us to check whether the univer- 
sal jump predicted in the KT picture is valid in systems 
defined by Eqn. ([TJ too or not. The present paper also 
explores the fact how disorder influences the properties 
of phase transition in these 2D systems. The effect of 
disorder on the KT transition has become relevant since 
the experimental observation of superconductor- insulator 
transition in thin disordered films [25l |26| . 

For the purpose of investigation, I choose the model 
defined by the interaction potential given by Eqn. ((T|). I 
have found that the transition is associated with a non- 
universal drop in the helicity modulus. I also find that 
when the number density of topological defects increases 
rapidly, disorder is introduced into the system. The ad- 
ditional parameter p 2 plays the role of disorder here. 

The variation of average defect pair density (p) with 
the dimensionless temperature T for lattice size L = 64 is 
shown in Fig. [T] The coupling constant J (in Eqn. ([1])) 
has been conventionally set to unity. The method for cal- 
culating average defect pair density and the simulation 
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FIG. 1. (Color online) Average defect pair density p plotted 
against dimensionless temperature T for L = 64 for various 
values of p 2 

techniques are discussed in Ref. 18]. A sharp variation 
of p as T increases through the transition temperature 
T c (p 2 ) is observed, p is found to show a sharp jump at 



T c (p 2 ), particularly for large values of p 2 . So, for strong 
enough nonlincarity, there is a sudden proliferation in the 
average defect pair density and the system under investi- 
gation behaves like a dense defect system for large values 
of p 2 . 

Next I plot the average defect pair density (p) as a 
function of the parameter p 2 , shown in Fig. O The 
plot is for three different system sizes at a temperature 
T = 1.1200, which is above the transition temperature 




FIG. 2. (Color online) Average defect pair density p plotted 
as a function of p 2 at T = 1.28 for three system sizes. The 
best fit corresponds to L — 64. The error bars are smaller 
than the dimension of the symbols used for plotting. 

of the model for p 2 = 50. The data for p versus p 2 are 
nicely fitted by the following expression 



p{T) = /9 max - a(T)exp(-£^/p 2 



(2) 



Eqn. ([2]) takes into account both vortices and anti- 
vortices, p increases monotonically with p 2 and saturates 
to a maximum value /? max f=a 0.28, which is independent 
of the temperature. The values of p max for the three sys- 
tem sizes are listed in Table HI There is no significant 
system size dependence of the parameters as is evident 



TABLE I. parameters for the fit of p(T) = p max 
a(T)exp(— 8y/p 2 ) for different L 



L 


Pmax 


a(T) 


8 


32 


0.2866 ± 0.002 


0.411 ±0.005 


0.301 ± 0.008 


48 


0.2876 ± 0.002 


0.406 ± 0.006 


0.297 ± 0.009 


64 


0.2878 ± 0.002 


0.406 ± 0.006 


0.296 ± 0.009 



from Table |TJ The maximum defect density (p max ) is usu- 
ally achieved in the high-temperature limit (T — ¥ oo) but 
here it is achieved in the high-p 2 limit (p 2 — > oo). There- 
fore it seems reasonable to interpret the parameter p 2 to 
play the role of disorder. In the high-p 2 limit, the system 
contains only vortex excitations. This means that in the 
high-p 2 limit, the system must be disordered even at very 
small temperature and consequently the transition tem- 
perature decreases with increase in p 2 which we observe 
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in Fig. [TJ In other words, as we increase the nonlinear- 
ity p 2 , the influence of disorder becomes stronger and a 
tendency of a first order transition with a sudden prolif- 
eration of topological defects develops. 

I now concentrate in the behavior of spin-stiffness (he- 
licity modulus). Helicity modulus, introduced by Fisher, 
Barber and Jasnow [27| . is a thermodynamic function 
which measures the "rigidity" of an isotropic system un- 
der an imposed phase twist. The free energy difference 
between the twisted and the periodic boundary condi- 
tions is proportional to the helicity modulus (7). Thus 
the general definition of helicity modulus may be re- 
garded as m 
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lim 2L 2 - 



F(w) - F(0) 



(3) 



where lo is the angle of twist and d is the spatial dimen- 
sion. In the present case, the twisted boundary condition 
being anti-periodic and d being 2, the definition of 7 sim- 
plifies to 



7 = 2 



F{n)-F(p) 



(4) 



Anti-periodic boundary condition is imposed only along 
one direction, say along X-dircction. Our calculation of 
7 involves a direct simulation of Eqn. (j4]) using multiple 
reweighting histogram method, a sophisticated MC tech- 
nique first proposed by Ferrenberg and Swendsen [29L f30| | . 
In the simulations, 10 7 MC steps per site were used for 
computing the raw histograms and 10 6 MC steps per site 
were taken for equilibration. 

The variation of 7 with temperature for various lat- 
tice sizes is displayed in Fig. [3] The transition is sig- 
naled by an abrupt decrease of the helicity modulus in 
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FIG. 3. (Color online) Helicity modulus 7 against tempera- 
ture (T) for various lattice sizes for p 2 = 50 with the errorbars 
shown (for three lattice sizes). 

the vicinity of the transition temperature as the tem- 
perature is increased and the drop at the transition gets 
steeper as the system size(L) is increased. It is also man- 
ifested that instead of the subtle and the smooth KT 



transition, the transition coincides with a non-universal 
jump in 7. The physical picture can be explained as 
follows. As p 2 increases, short-scale fluctuations (rota- 
tion of separate spins by large angles) are favored over 
long-wavelength fluctuations (spin waves and vortices) 
and when this happens disordering is induced at a lower 
temperature than the KT transition temperature, but 
the vortex-vortex interaction still remains stronger at 
that lower temperature, thus making the helicity mod- 
ulus jump non- universal. The ratio of 7/T at T = T c 
for different L is listed in Table [TTJ I point out that 
Minnhagen Q also showed the possibility of a KT tran- 
sition in a 2D CG with a non-universal jump in 7. 



TABLE II. 7/T values at T = T c for different L 



L 


16 


32 


48 


64 


80 


96 


7/T 


0.5742 


0.6049 


0.6569 


0.7192 


0.6971 


0.6454 



I have also computed the temperature derivative (t) of 
the helicity modulus. The internal energy difference un- 
der anti-periodic and periodic boundary conditions gives 
the derivative of helicity modulus in the form 



1 d \R (RW ^ a ~ 

T= 2dp [Pim = Ifi 



(E) p 



(5) 



The MC data for the right hand size of Eqn. ([5]) as 
a function of temperature for different lattice sizes are 
shown in Fig. |4j The transition is manifested by a huge 
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FIG. 4. (Color online) Derivative of |/?7(/3) as a function of 
temperature T for different lattice sizes for p 2 — 50 with the 
errorbars shown (for three lattice sizes). 

peak height in r and the data display a divergent behav- 
ior with increasing L, indicative of a discontinuous jump 
in 7 in an infinite lattice. 

Now I present the FSS of r. Since r is a response 
function like specific heat (C v ) or susceptibility (x), it is 
expected to show an identical behavior in scaling as for 
C v or x- From Fig. [51 where the maxima of r are plotted 
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against L 2 , it is clear that the peak heights of r scale as 
L d which confirms the first order nature of the present 
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FIG. 5. (Color online) The peak heights of r plotted against 
L 2 with the linear fit represented by the straight line. The 
error bars for most points are smaller than the dimension of 
the symbols used for plotting. 

model for p 2 — 50. We recall here that for first order 
transition standard scaling rule for C v goes like C v ~ L d 
31|. 

In this Communication I have used extensive MC sim- 
ulation to show that for strong enough nonlinearity (i.e., 
for large values of p 2 ) in the interaction potential of 
Eqn. |T]), there is a sudden proliferation of topologi- 
cal defects which makes the system disordered. Conse- 
quently the transition is associated with a discontinuous 
non-universal jump in the helicity modulus. Thus our 



simulation has given some support to the idea that the 
type of phase transition in thin superconducting films 
may be changed due to influence of disorder. As high- 
T c superconducting films are believed to have a irregular 
structure, it seems reasonable to relate the non-universal 
jump to disorder. 

However, some studies [lfllHUIHil; mostly based on RG 
analysis of Migdal-Kadanoff type contested the first order 
nature of the transition in the model defined by Eqn. ((T|) . 
Since renormalization arguments hold good only for small 
disorder, the possibility that disorder may change the 
nature of phase transition always remains there. Perhaps 
this is the reason behind the different interpretation of 
results by the authors of Ref . [lg, [32J, |33J . 

Finally, the present work could shed light on the nature 
of 2D melting which remains controversial for decades. 
Experimental works and numerical simulations favor a 
KT-like transition in some cases and a discontinuous one 
in others j34|. After all, it is possible that the nature 
of the melting transition in 2D depends on the specific 
system and the parameters of the model which in turn 
translate into different values of the nonlinearity param- 
eter p 2 . 

I end this paper with a comment. It is observed in 
Fig. [3] that the graphs for helicity modulus (7) against 
temperature for different lattice sizes intersect at a point 
which is the transition temperature of the model (within 
an error of 0.03%). I offer no explanation for this inter- 
esting result and this issue is left for future research. 

The author is grateful to S. K. Roy for many fruitful 
discussions and a number of suggestions after critically 
reading this manuscript. 
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